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Abstract. - We show the existence of three well defined time scales in the dynamics of 
wormlike micelles after a step between two shear rates on the stress plateau. These time scales 
are compatible with the presence of a structured interface between bands of different viscosities 
and correspond to the isotropic band destabilization during the stress overshoot, reconstruction 
of the interface after the overshoot and travel of a fully formed interface. The last stage can 
be used to estimate a stress diffusion coefficient. 



Introduction. - Depending on the type and concentration of surfactant molecules and 
added salt, solutions of surfactant wormlike micelles have shear thinning or thickening be- 
havior under shear flow. Unlike most fluids, wormlike micelles often have non-analytic flow 
curves with sharply-selected plateaus along which strain rate or stress may change discontin- 
uously. In the well documented case of shear thinning solutions the usual explanation of the 
constant stress plateau is shear banding [^H,^, i.e. a separation of the material into bands 
of different viscosities, triggered by a constitutive instability (such as an isotropic-to-nematic 
transition |^). As shown recently |§, ||, |lO|, |l^, |l^ , the stress selection and history indepen- 
dence of shear banding can be explained using the inhomogeneities of the relevant mesoscopic 
order parameter (polymer stress), i.e. by incorporating "diffusive" terms in the constitutive 
equations. Order parameter diffusion was introduced long time ago by van der Waals in 
^ I the so-called "gradient theory" of the gas- liquid interface [Q , and is obligatory in phase field 
• ' models for pattern formation. Notwithstanding a few attempts to deal with inhomogeneous 
rN " stresses the same concept has not obtained full acceptance in the rheological community. 

. While one might argue that diffusion terms are negligibly small, these non-perturhative terms 
■ ■ ■ ' resolve stress selection even for infinitesimal values However, a small diffusion coeffi- 

cient should also imply a slow approach to steady state; the main purpose of this letter is to 
demonstrate these long time scales experimentally and relate them to simple model diffusive 
behavior. 
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Table I - Surfactant systems used in this study; the surfactant was 0.3 M CTAB. r and Go, are the 
Maxwell relaxation time and plateau modulus, ji^Sj are the start and width of the constant stress 
plateau, (ri)i=i,3 are the three time scales discussed in the text, and a* is the plateau stress. 

No Salt T [°C] Go [Pa] r [s] n [s] ra [s] rg [s] rj-/ r^7 a* jGo 

1 1.79M NaNOs 30 232 0.17 0.2 2.4 30.6 0.85 19 0.64 

2 0.405M NaNOs 30 238 0.17 0.2 1.8 25 1.27 19 0.66 

3 0.3M KBr 34 235 0.16 0.2 3.7 9.3 1.12 80 0.66 



Shear banding involves spatial inhomogeneity and several temporal stages. Light polariza- 
tion probes the local micellar orientation, while rheology detects the molecular stress. In this 
work we shall calculate the stress transients using a theoretical model. This will be compared 
to birefringence measurements that are assumed to probe the state of molecular orientation 
and hence stress (but see Ref. [|3|) Rather than the typical start-up transient experiment, we 
consider the simpler experiment of a step between two fixed values of the shear rate in the 
banded regime. Small steps should induce less drastic changes in the fluid while still remaining 
in the non-linear regime, and hopefully yield more controllable results. The transient features 
will be shown to be intimately related to the dynamics of the interface between the bands. 

Experiments. - The surfactant solutions used are summarized in Table ||. The stress 
response was measured using an RFS Rheometrics Scientific controlled strain rate rheometer 
in Couette (radii 24.5,25 mm) and cone-plate geometries. Linear response is of the Maxwell 
type with almost identical relaxation times t and moduli Gq for the three solutions. In the 
non-linear regime the coexistence plateau width Sj for solution 3 is larger than for solutions 
1 and 2, while solutions 1 and 2 have a sloped stress plateau that roughly follows a power 
law tJ ~ 7", with a < 0.1; a slope could indicate concentration differences between coexisting 
states . Upon a step increase of the shear rate between banded states the stress increases 
to a maximum, and then decreases monotonically (or sometimes by a small undershoot and 
a monotonic increase). Three relaxation times (ri w r w 10^^ s, T2 ~ lOr wis and T3 « 
IOOt w 10s) follow successively after the overshoot until steady flow is reached (Fig. ^). The 




Fig. 1 - Stages and characteristic times (in r units): a) total stress during the step 7 = 10 — * 20 s 
for CTAB/KBr; b) simulated total stress during the step 7 = 1.4/r 2.8/r using the d-JS model; 
c) Extinction angle during the step 10 20 s"^ for CTAB/KBr. 
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three time scales are well separated and the result is reproducible for all solutions and both 
cone-plate and Couette geometries. We filmed the step shear rate experiment for the third 
solution in a Couette geometry with a slightly larger gap (radii 24, 25 mm) between crossed 
polarizers to extract the average extinction angle x- The kinetics of x resemble that of the 
shear stress, displaying time scales similar to the second and the third time scales found in 
the rheology (Fig. |l|c). The first, unresolvable, time scale is shorter than the interval between 
successive video frames. The second time scale is well resolved, although shorter by a factor 
of two than the equivalent time in the stress rheology measurements. The third time scale is 
buried in noise and regression does not provide significant results. 

Theory. - The momentum balance is p {dt + V) v = V-T, where p is the fluid density 
and V is the velocity field. The stress tensor T is given by T = —pi + 2riD + S, where the 
pressure p is determined by incompressibility (V-t; ~ 0), 77 is the "solvent" viscosity, S is the 
"polymer" stress, and D is the symmetric part of the velocity gradient tensor (^v)ai3 = daVp. 
The non-Newtonian "polymer" viscoelastic stress S is assumed to obey the diffusive Johnson- 
Segalman (d-JS) model 

{dt + v-V)ll~ (rJS - E f2) - a (DS + HD) = PV^S + 2/iZ)/r - S/t, (1) 

where $7 is the anti-symmetric part of Vw, p, = GqT is the "polymer" viscosity, r is a relaxation 
time, and T) is the diffusion coefficient. The "slip parameter" a (describing the non- affinity 
of the deformation) is necessary to reproduce a non-monotonic constitutive curve, and the 
added diffusion term was shown to resolve stress selection ^ . 

The initial dynamics is governed by inertia; within a very short time tm ~ pL^ hi [tm = 
10~^ s for rj/ p — 0.01 m^s~^ and a gap L = 1 mm) momentum diffuses across the gap and the 
momentum balance becomes T = const. The subsequent slower dynamics is controlled by 
the viscoelastic response of the fluid. In a planar geometry with v — v{y)x, our constitutive 
model leads to a system of reaction-diffusion equation^Tl: 

where 7 is the shear rate, S = Y^xy, and is a combination of the polymer normal stresses, 
Tixx and Syy. S^W are the order parameters of the transition [S is small in the nematic 
(N) band, and large in the isotropic (I) band). They can diffuse across stream lines with 
diffusion coefficient I?, and relax in the linear regime within the linear (Maxwell) time r. The 
non-linear reaction terms Cs = i{Go ~ 1^) and Cvy = 7S' can be straightforwardly derived 
from Eq. (0) [||lO|. 

The local momentum balance for the shear stress a = T^y is 

o = S + eGoT^. (3) 

where e — rj/ p. The dynamics of Eq. ^can be understood with the aid of two (local) dynamical 
systems: 

S^-S/T + Cs{i.S,W), W = -W/r + Cw{i,S,W)- (4a) 

S^-S/T + Cs{{(J~S)leGoT,S,W), W = -W/t + Cw{{(t - S)/eGQT,S,W) (4b) 

(^)Any non-monotonic difTcrential constitutive model (e.g. Gates Jl?)], Doi-Edwards [|l6|,|l5|]) with diffusion 
terms leads to a similar equation set. 
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Fig. 2 - a) Relaxation times to I and N attractors at constant stress, tI ,t'^. b) Thick line: theoretical 
flow curve (the negative slope,unstable branch should be replaced by the constant stress plateau). 
Thin lines: trajectories of the coexisting bands near the walls, {si}, i = f , 3 are the starting points of 
the three stages discussed in the paper. 



where S = dS/dt. System (Q) describes the dynamics along a streamline at prescribed shear 
rate; in this case a changes proportionately to S according to Eq. (||) . System (^) describes 
the dynamics along a streamhne at constant total stress a. 

The two dynamical systems have the same fixed points (since, for homogeneous steady 
flow, 7 and a are related by Eq. ^): stable fixed points (attractors) representing the bands 
I and N, and an intermediate unstable saddle fixed point. Coexistence of bands at common 
total stress is possible only for a £ [ai{e),a2{e)\. Linearizing systems ( pa| ) and (^) about 
the fixed points yields the dominant relaxation times of the attractors, t/ and tjv- These 
are different for the two dynamical systems, denoted at constant shear rate by t] ^ and at 
constant stress by r/jy. For the JS model t] = t]^ = t for all 7. is close to r, while 
rf is larger than and diverges as cr — > a2(e) (Fig. ^). This divergence is consistent with 
Ref. 1^: controlled stress experiments have much longer relaxation times on the metastable 
extension of the high viscosity branch above the constant stress plateau. 

Consider an initial banded steady state, with average shear rate (7) . Suddenly increasing 
the average shear rate to 72 > (7) produces a stress overshoot because the amount of high 
viscosity I band is too large; both I and N bands then become unstable, and the stress can 
decrease by producing more of the low viscosity N material. Numerical simulation (Figs.|l|, ^ 
shows that this occurs in three stages: 

1. Band I destabilization — During this stage the I band tries a direct passage toward the 
nematic band N, Fig.^). Complete transformation is forbidden by the average shear rate 
constraint and the 1 band stops before reaching the basin of attraction of the steady N band. 
A representation of the subsequent kinetics in the (7, a) plane (Fig. |^) shows that the N band 
almost follows the steady flow curve, while the I band evolves at constant shear rate. Thus, the 
characteristic time ti « « r is controlled by the 1 band dynamics. The total shear stress at 
the end of this stage depends on the final position of the interface and is sensitive to the details 
of the constitutive model. If this value is below the plateau one finds an undershoot (as in 
the numerical simulation Fig. |l|b) , and otherwise the subsequent stress evolution is monotonic 
(Fig. |l^). Although we did not succeed in reproducing monotonic evolution using the JS 
model, as long the I band evolves at a constant shear rate, the magnitude of the characteristic 
time should not be affected by the presence or absence of an undershoot. 

2. Interface reconstruction — At the end of stage 1 the interface separates an unsteady I 
band close to the unstable saddle point from a nearly stable N band. The part of the profile 
closer to the I attractor will evolve toward this one, while the other part approaches the N 
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Fig. 3 - Simulation of the order parameter profiles for the three stages (s: start, e: end): a) destabi- 
lization (0 < t < ti); b) reconstruction {t\ < t < 12); c) travel {t2 <t < 00). 

attractor. This reconstructs the interface in a more advanced position, stabihzes the bands, 
and increases the contrast between them. Interestingly, there is a spatial position at which 
S and W practically remain fixed at their saddle fixed point values. During this stage (7) 
is constant (because it is imposed) and {S) is almost constant (because of the compensating 
evolutions of the two bands), so according to Eq. (||) the total stress variation is small, and 
the characteristic time (controlled by the I band) is T2 « rf . This time exceeds the linear 
viscoelastic time r (see Fig. and depends on how close the stress at the end of stage 
1 is to the spinodal limit cr2(e), and on the quantitative details of the curve in Fig. |, aU 
of which are sensitive to the constitutive model. The analysis suggested by the numerical 
experiment is confirmed by the birefringence measurements. The sequence of images in Fig. ^ 
show the gap of the Couette cell filmed between crossed polarizers during stage 2. Although 
we can not quantitatively compare Figs, ^b) and ^ (the relation between the transmitted 
intensity and the order parameter is unknown and sure to be non-linear), the sharpening of 
the contrast corresponding to the interface reconstruction is visible. The difference between 
the characteristic times for the extinction angle and rheology (Figj^) could be due to the 
different Couette cell gap widths. 

3. Interface travel — The instability and reconstruction of the interface in the first two 
stages is ensured by the reaction terms of the Eq. (^ , ending when a sharp interface between 
stable bands is fully formed. This interface could have a non-zero velocity if it forms in 
at a position corresponding to a stress value above or below the plateau stress a* . "Front 
propagation" over the small distance toward the final equilibrium position is then controlled 
by T) (this distance is too small to observe by birefringence). Because of the undershoot in 
the numerical simulation the sign of the displacement during stage 3 is opposite to the one in 
the first stages (Fig.^). The characteristic time T3 for this stage follows from the velocity c 
of the sharp interface close to steady state which is independent of the presence or absence of 
the undershoot |0 (see below). 

Let us consider a single sharp interface, at a position r inside the gap. At imposed shear 
rate, the lever rule (7) = j^i'^ {(y) + (1 — ■j;)Y{'y) relates a and r and leads to 
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Fig. 4 - CTAB/KBR: 7 = 10 ^ 30s~^ jump: birefringence images and profiles (averaged between 
the two vertical lines on the film) corresponding to the second time scale, showing the interface 
reconstruction. The moving (inner) cylinder is at r = 1.0 and the fixed (outer) cylinder is at r = 0.0. 



where 5'^ = jN—ji is the width of the plateau, and rjj = da/d^\^^. We consider (7)— 7/ 1^7, 
(57 S> 7/ (as in the experiments) so that r/^Y ~ rj^Sj (true for piecewise linear flow curves 
and obeyed well by the JS model) . 

We showed previously that the velocity c of the interface is a function of the total shear 
stress tJ only, and that c — when a = a * |,0 which via the lever rule corresponds to a 
unique stable interface position r*. Thus, close to this position the equation of motion of 
the interface is: 

i-.i.,^t{^)jr-., (e) 

Using Eqs. (||,^, and the derivative ^\^_^, = KGoy^, we find the solution r — r* = 
(r(0) — r*)e~*/'^^, with characteristic time 

L 1 + [(7) ~7/]/7J 

where the dimensionless constant K depends on the particular constitutive model. 

Eq. (|^) implies that a fully formed interface equilibrates faster in systems with larger 
plateaus ^7, such as CTAB/KBr. In such cases, Eq. (|^) implies larger stress variations and 
thus larger interface accelerations for the same position variation. This is compatible with the 
shorter ra in Table |. For simplicity, Eq. (||) was for a planar geometry; in cylindrical Couette 
flow a slight correction (negligible for the thin gaps we consider) leads to a smaller t^. 

Using the experimental value of and Eq. |^, we can estimate T). In order to do this 
we need the value of K. In the JS model, while Go and t are measurable, the two free 
parameters e and a determine the function K. Nevertheless, tGqS'^ /(J* = /i(e), ji/Sj = /2(e) 
and K/{t5'^) = /3(e) are functions of e only, given to a good approximation (> 80%) by 



/i(e) = i- / /2(e) = 3| ^ |eVi(e)- From either /i, or ^ and 

experimental data, one can estimate e (the average of the two values is given in Table || 
and then /a gives K (Tabic Thus we do not need the value of a. Microscopically, we 
expect T) = C,"^ /t where C is the stress correlation length. In dilute solutions this should be 
the micelle gyration radius H, while in concentrated solutions a reasonable candidate is the 
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Table II - Stress diffusion estimates using the JS model. D ts obtained from the values of T3 (table 
and Eq. VA with e, K estimated using 7/ and 5^ . 



Sample 


L [mm] 


(7 - ii)hi 


e 


K 


D [m^s-i] 


C [nm] 


i [nm] 


1 


0.3 


4 


0.023 


0.10 


11. 10"'^ 


44 


26 


2 


0.3 


2.33 


0.028 


0.12 


11. 10"^^ 


43 


26 


3 


0.3 


1.85 


0.007 


0.14 


3. lO"^'^ 


22 


26 



mesh size which can be estimated from Go ~ kT/£_^. The results are presented in Table 
The stress correlation length is of order the mesh size, which is reasonable; however there is 
still no theory for such a diffusive term in concentrated solutions. 

To conclude, a general dynamical systems analysis of the d-JS model provides plausible 
explanations for the observed time scales, and consistent estimates of the stress diffusion 
coefficient V. Nevertheless, neither the d-JS model, nor reptation-retraction-reaction models 
[p^,[r^] can provide perfect fits of transient stress curves in the nonlinear regime. Also, it 
it conceivable that concentration differences between the bands could influence the observed 
time scales. The difference of the values of V between samples 1,2 and 3 (Table ||) could be 
a concentration effect, consistent with the different slopes of the flow curves plateaus (tilted 
for 1,2, almost horizontal for 3). 
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